Geometric Mechanics of Curved Crease Origami 
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Folding a sheet ol paper along a curve can lead to structures seen in decorative art and utilitarian 
packing boxes. Here we present a theory for the simplest such structure: an annular circular strip 
that is folded along a central circular curve to form a three-dimensional buckled structure driven 
by geometrical frustration. We quantify this shape in terms of the radius of the circle, the dihedral 
angle of the fold and the mechanical properties of the sheet of paper and the fold itself. When the 
sheet is isometrically deformed everywhere except along the fold itself, stiff folds result in creases 
with constant curvature and oscillatory torsion. However, relatively softer folds inherit the broken 
symmetry of the buckled shape with oscillatory curvature and torsion. Our asymptotic analysis 
of the isometrically deformed state is corroborated by numerical simulations which allow us to 
generalize our analysis to study multiply folded structures. 



Over the last few decades, origami has evolved from an 
art form into a scientific discipline [jj, [2| , made possible 
in large part by methods for the mathematical analysis 
of folded surfaces. Though there is artistic and techno- 
logical precedent for folding paper along curves [SHjj, as 
seen in Japanese paper boxes, and the ubiquitous Mc- 
Donald's fries box, how these geometrical structures get 
their physical shape is poorly understood. In this letter, 
we consider the geometric mechanics of the simplest case 
of a curved fold: an annular elastic sheet folded along a 
closed, circular curve which buckles out of the plane to 
resolve a fundamental incompatibility between the folded 
geometry and the ensuing mechanical stresses. 

Qualitative experiments with a complete circular an- 
nulus of paper having a concentric, circular crease show 
that folding buckles the crease into a saddle (Fig. Q}t), 
while the same crease along a cut annulus remains planar 
(Fig. [Tp). This behavior is a consequence of a fundamen- 
tal incompatibility between the geometry of the fold and 
the stretching elasticity of the sheet. As we will see, and 
it is apparent in Fig. [lb, the sheet responds to folding by 
wrapping around itself to eliminate in-plane mechanical 
stresses. The closed annulus, on the other hand, can ex- 
pel these stresses by buckling. In the limit of where the 
thickness of the sheet is much smaller than the width, 
which is itself smaller than the length of the crease, the 
shape that arises is a balance between the bending energy 
of the sheets on either side of the crease, the energy at 
the crease itself, and the geometrical constraints arising 
from the sheet's closed topology. Topological constraints 
are also crucial to shape and mechanics in the small de- 
formations of shells and non-Euclidean plates Q . 

We consider an annulus of uniform thickness t, width 
2w folded along a central circular crease of radius r 
(t <C w < 2nr). In the deformed state, the crease is a 
space curve parametrized by arc length s, with curvature 
k(s) and torsion r(s), and the surfaces on either side of 
it come together at a finite fold (or dihedral) angle 0(s). 



Assuming isometric deformations away from the crease, 
the mid-surface of the sheet on either side of the crease is 
developable. Then any point on it can be characterized 
in terms of a set of coordinates (s,v), corresponding to 
the arc length and the generators of the developable on 
the inside and outside of the crease, g± (Fig. [It) , with 
the coordinates: X±(s, v) — X(s, 0) + vg±(s). For devel- 
opables, the generators must also satisfy the condition 
that g±(s) x dg±(s)/ds is perpendicular to the crease 
jlQ]. Since folding does not induce in-plane strains, the 
projection of the crease curvature onto the tangent plane 
on either side of the sheet must remain 1/r. This leads 
to two geometrical conditions [B| that relate the dihedral 
angle of the crease to its spatial curvature and the angle 
of the generators of the developable surface on either side 




FIG. 1: (color online) A photograph of the model built by 
cutting a flat annulus of width 2w from a flat sheet of paper 
with central circle of radius r. (a) Folding along its center 
line buckles the structure out-of-plane. However, if we cut 
the annulus, (b), the structure collapses to an overlapping 
planar state with curvature given by equation {TJ. (c) The 
inset shows a cross section of the fold, where the right and 
left planes, X+ and X_, define the dihedral angle 9. 
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where n/r and r/r are the curvature and torsion of the 
crease respectively, 7± is the angle between the unit tan- 
gent vector of the crease d~K(s)/ds and the generator. 
We see that k(s) > 1, with equality only when 9 = n. 
For a circular crease concentric with a circular annulus 
of constant dimcnsionlcss half-width u) = w/r, we find 



The equilibrium shape of the curved crease results from 
minimizing E = Ef, + E c and is characterized by three 
parameters: the scaled natural width of the ribbon uj, the 
natural dihedral angle between the two surfaces adjoining 
the crease 9o and the dimensionless crease-surface energy 
scale cr, subject to appropriate boundary conditions. For 
example, an open circular crease has free ends and thus 
prefers to remain planar with t = since non-planarity 
would increase both the curvature and torsion (l3| . A 
closed crease, however, is frustrated by geometry, forcing 
it to buckle, a fact that follows from the inequality k = 
1/ sin(0/2) > 1 when 9 < tt which requires J dt; k > 2ir, 
and is incompatible with a planar crease with r = 



^(0 = ±sin 7± (0 T y w 2 ^ 2w + sin^ 7± (0 (3) 

to be the dimensionless distance to the boundary along a 
generator leaving the crease from a point labeled by the 
dimensionless arc length £ = s/r. 

The energy of the sheet is the sum of the energy of 
deforming the sheet on either side of the crease and that 
of the fold that connects them. Since the creased folded 
surface is piecewise developable, the energy per unit sur- 
face is proportional to the square of the mean curvature 
111. The mean curvature on either side of the sheet is 
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where (.)' = d/d£(.). Then the energy of each surface 
Eb = B J Q 27r f?* H±dvd£,, where B is the bending stiff- 
ness of the material of the sheet. Carrying out the inte- 
gral along the generators, v, explicitly leads to the fol- 
lowing scaled bending energy for the two surfaces 
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We see that © is determined entirely in terms of the 
geometry of the crease. To model the fold itself, we use 
a phenomenological energy functional measuring the de- 
viation of from an equilibrium angle 9g, which we 
assume to be constant, so that the scaled crease energy 
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where a = Kr/B is the ratio of the crease stiffness K 
and the bending stiffness B. This energy reduces to a 
simple quadratic expression in the difference 9 — 9o when 
9 ~ 6*o ; although the precise form of this term does not 
affect our analytic results, it conforms to our numerical 
model 0. 
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FIG. 2: (a) Perturbative fold of width u = 0.1 and a = 2/y/3 
shaded by mean curvature. The generators are indicated by 
the lines on the surface. The inset shows the dimensionless 
torsion and curvature of the crease, (b) A simulated fold of 
width uj w .0994 shaded by local area change relative to the 
flat state. 

Though geometrical constraints induce buckling, the 
resulting fold shapes are determined by minimizing the 
total elastic energy consisting of contributions from the 
sheet ([5]) and the fold ([5]), expressed entirely in terms 
of the curvature and torsion of the crease [1J, [l5 1 . For 
relatively narrow, but stiff, folds i.e. w«l and a 1 
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that are weakly folded, i.e. so that the dihedral angle 
9 ~ 7T, and thence e = l/sin(0 o /2) - 1 < 1. Then, 
we find that the total scaled energy E = (E b + E c )/B 
simplifies to [l3[ 

E*l" «{£(„- 1-. + (7) 

in terms of the scaled curvature k and torsion r. We 
see that as a — > oo, the rescaled curvature k — ¥ 
l/sin(^o/2) = 1 + e, the prescribed curvature. The 
minimal energy crease shape, therefore, minimizes r 2 
subject to the constraints of fixed length and curva- 
ture. In this limit, the Euler-Lagrange equations become 
[t" + (1 + e) 2 r]' w at constant curvature [13]. If e = 0, 
corresponding to a dihedral angle 9q = ft the solution to 
these equation is infinitely smooth. Otherwise, a solution 
of continuity class C may be obtained to these equations 
with k = 1 + e and oscillating torsion, 
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The absolute magnitude of the torsion tq is then de- 
termined by the condition that the curved fold has arc 
length 2irr, and consistent with the four-vertex theorem 
for closed convex space curves, there are four points with 
vanishing torsion 16 1. 

To go further, we can carry out an asymptotic analysis 
of the Euler-Lagrange equations obtained by minimizing 
E = Eb + E c must be performed by expanding the shape 
of the crease around a planar curve of constant curvature, 
kq. Following 3, 15 1, we write k = kq + 5 k and r = St 
and compute the Euler-Lagrange equations. To lowest 
order, we obtain an algebraic expression determining the 
ideal curvature of the crease, kq for arbitrary a, e and 
u> [IH, the curvature of an incomplete or severed planar 
annular fold with zero torsion. To next order, we find 
that both the curvature and torsion oscillate; a typical 
analytical solution is shown in Fig. [5^, with the inset 
showing the oscillating torsion vanishing at the extrema 
of curvature (see Fig. [IJ. Here the overall amplitude 
of t is chosen to close the curve, with 9(0) and 9(n/2) 
parametrizing the solutions [l3| . 

These qualitative features are also confirmed by di- 
rect numerical minimization of the energy of a triangular 
mesh model for the curved origami structure in which 
each edge is treated as a linear spring and adjacent tri- 
angles in each sheet have collinear normals and , with the 
scaled ratio of the bending stiffness to the stretching stiff- 
ness B/Sl 2 pa 10~ 3 . However, adjacent triangles across 
the crease prefer a fixed, non- planar dihedral angle [l2j ]. 
The presence of a small but finite extensibility of this 
model implies that our simulations relax the isometry of 
the folding process and thus allow us to capture how ex- 
tension and shear arise in wide folds (Fig. [2]b) . We find 



that the extensional and shear strains typically localize 
where the mean curvature becomes large, consistent with 
our isometric analytic theory (shown in Fig. [2^). 
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FIG. 3: (a) Angle differences \6{n/2) - 0(0) | as a function 
of uj with 9q — 2tv/3. The red curve (diamonds) are com- 
puted from first-order perturbation theory with a — 2/\/3 
and 6o = 2-7T/3. The numerical simulations (blue - dashed 
lines) are shown with a — 2\/3 (circles) and a = 160/\/3 
(squares) and are compared with a nonperturbative varia- 
tional ansatz (green - solid line), Km and Tp), described in 
the text with a — x/3/40 (circles) and a — V3/2 (squares). 
Corresponding energy landscapes, as a function of 8(0) and 
8(n/2) respectively, are shown for (b) w = 0.01, (c) 0.05, and 
(d) 0.1, with energy minima drawn as white dots. 



Moving beyond the simple asymptotic theory for nar- 
row folds, we consider the dependence of the solution on 
the scaled width by using the perturbative shapes as a 
variational ansatz in the exact expression for the energy 
Eb + E c . Since the shapes have a 4- fold symmetry, we 
expect to see a coincidence between 9(0) and 9(n/2). In 
Fig. [3h., we plot |0(O) — 9(n/2)\ for the minimal energy 
configuration as a function of the scaled width w. We 
see that when uj < 0.1, annuli with large a have a nearly 
constant dihedral angle around the entire length of the 
fold, with 9(0) pa 9(ir/2) for the narrowest fold widths. 
However, for small er, the energy minimum generically 
has 9(0) 9(ir/2); this discrepancy increases with the 
scaled width ui. Plotting the corresponding energy land- 
scape in Fig. [3p-d for some representative values of u, 
we see that the energy contours develop forks because a 
range of 9(0) and 9(n/2) are forbidden by the geometric 
constraints that the generators of our two surfaces can in- 
tersect only outside the actual surface, else the bending 
energy diverges. 
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To avoid the intersection of generators inside the outer 
surface requires 



, sin 7+ 
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which expression reduces to | r' | < (1 — uj) cot (6/2) /u>, at 
points where r = 0. Similarly, to avoid the intersection 
of the generators on the inner surface inside the inner 
boundary requires the discriminant in equation ([3]) to be 
positive, implying a bound on the torsion, 
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These geometrical bounds restrict the range of allowed 
torsion and thus the buckling of the crease. As a con- 
sequence, wide folds will become resistant to deforma- 
tions as the sheet quickly reaches a regime where the 
generators start to nearly intersect in the neighborhood 
of £ = 7r/2. In Fig. [3ji, this is manifested by the presence 
of large forks carved out by the forbidden configurations. 
Since the energy minima occur close to the singularities, 
our perturbative expansion of the shape is approximate 
at best. However, even at intermediate widths, where 
the perturbative expansion should be at least qualita- 
tively valid, the bifurcation of the minima show up in 
the shadows of the prominent forks observed in Fig. [3J1. 

These calculations suggest a second improved ansatz: 
«(2) = K o + K\ cos(2£) and T( 2 ) = r [sin(2£) + ^sin(4£)], 
choosing To to close the fold and f] to minimize the en- 
ergy. When r] = we now find very good agreement with 
the perturbative ansatz previously considered. However, 
we find that rj w —0.45 for large widths, which lowers the 
maximum of the torsion and better satisfies the singular- 
ity bounds in Eqs. (|9l- IT0| . Using a as a fitting parame- 
ter, we see that 0(ir/2) — 0(0) agrees quite well with the 
numerical solutions for small uj and only diverges from 
numerical simulations for large widths, around U) ~ 0.08 
as shown in Fig. [3^. 

Finally, we consider structures built from multiple, 
concentric folds (Fig. SJ. Again, the large penalty for 
stretching leads to strong geometric constraints connect- 
ing the curvature and torsion of the crease to the dihedral 
angle of the fold given by equation ([1]), leading to self- 
similar creases and folds. The generator on the outside 
of a crease must coincide with the generator on the in- 
side of the next crease, fixing the angle 7+ on the inside 
of the next crease, while the torsion determines the rela- 
tionship between 7+ and 7_ , so that the direction of the 
next set of generators emerges. This procedure follows 
from the first crease to the last crease until we reach a 
boundary or the generators cross. Our numerical simula- 
tions confirm this and further show that multiple creases 
have little stretching content and do buckle rigidly. 
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FIG. 4: (color online) (a) A plastic model of six circular folds 
generated by perforating the folds at equal intervals with a 
laser cutter. The ratio of outer to inner boundary is 2. (b) 
A simulation result with the same planar geometry as the 
plastic model in shaded by local area change. The multiple 
folds simulation has the same magnitude order of area change 
as that of the single fold simulation. 



Our study of curved crease origami shows that a con- 
sequence of the fundamental frustration between fold- 
ing along a curve and the avoidance of singularities and 
in-plane stretching imposes geometric constraints on the 
shape that are reflected in a bifurcation of the curvature 
of a closed crease of large width. Indeed, the coupling be- 
tween shape and in-plane stretching endows these struc- 
tures with a stiffness and response that is unusual, as we 
have demonstrated in the simplest of situations - a closed 
circular fold. Moving forward, our approach may be gen- 
eralized to more complex curves with variable dihedral 
angles in folded structures with curved creases and thus 
sets the stage for the analysis and design of these objects. 
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